Magnetic fields and the outer rotation curve of M31 
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ABSTRACT 

Recent observations of the rotation curve of M31 show a rise of the outer part that can not 
be understood in terms of standard dark matter models or perturbations of the galactic disc 
by M31's satellites. Here, we propose an explanation of this dynamical feature based on the 
influence of the magnetic field within the thin disc. We have considered standard mass models 
for the luminous mass distribution, a NFW model to describe the dark halo, and we have added 
up the contribution to the rotation curve of a magnetic field in the disc, which is described by 
an axisymmetric pattern. Our conclusion is that a significant improvement of the fit in the outer 
part is obtained when magnetic effects are considered. The best-fit solution requires an amplitude 
of ~ 4 /iG with a weak radial dependence between 10 and 38 kpc. 

Subject headings: galaxies: individual (M31) , galaxies: magnetic fields 



1. Introduction 

Recent high sensi tivity measurements of the ro- 
tation curve of M31 ( Chemin et al ■l2009HCorbelh et al 
|2010[ hereafter C09 and CIO respectively) suggest 
that it highly rises at the outermost part of the 
disc of M31 (r > 25 - 30 kpc) (see Figure p. 
This behaviour cannot be considered an exception. 
Similar outer ris ing rotation curves can be f ound in 



other galaxies. iNoordermeer et al.l ([20071 ) stated 



that "in some cases, such as UGC 2953, UGC 3993 
or UGC 11670 there are indications that the rota- 
tion curves start to rise again at the outer edges 
of the HI discs", suggesting follow-up observa- 
tions of higher sensitivity to investigate this fact 
in more detail. Other examples co uld be found 



i n som e rotation curves provided by iSpano et al 
(|2008[ ). in particular the curves of UGC 6537 
(SAB(r)c), UGC 7699 (SBcd/LSB), UGC 11707 
(Sadm/LSB) or UGC 11914 (SA(r)ab) measured 



by means of Fabry-Perot spectroscopy. Other 
poten tial examples could be NGC 1832 and NGC 
2841 ( Kassin et a,l.ll2006h . even with a larger error 



and ESO 576-G51 and ESO 583-G7 ( Seiear et al 



20051) . even if the outer rising is small. Potential 



nearby candida tes could be NGC 31 98, DDO 154 

~ 20081) . Therefore, 



or NGC 7331 (d e Blok et al 



outer rising rotation curves are not uncommon 
in spirals, being the case of M31 the most repre- 
sentative. This puzzling dynamic feature requires 
a theoretical explanation. A detailed study by 
C09 showed that the NFW, Einasto or pseudo- 
isothermal dark matter halos fail to reproduce the 
exact shape of the rotation curve of M31 in the 
outer region and found no differences between the 
various halo shapes. Moreover, they found new HI 
structures as an external arm and thin HI spurs 
in the outskirts of t he disc. These spu rs have 
been also observed by iBraun et al. ( 20091) . Thus, 
a primarily explanation of this gas perturbations 
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could be the i nteractions with the M31's satellites 
asNGC 205 (iGeehan et all 120061: ICorbelli et al 
The main problem of matching the prop- 



erties of the giant stellar strea m observed in the 
south of M31 (jlbata et allbOOll ) is that the orbital 
of the companion that pro duce the stream is not 
constrained satisfactorily ( Fardal et al. 20061 ) or 



high velocities for the r adial orbits are found (see 
e.g. lHowlev et alJliooih . Although the northeast- 
ern (NE) and southwestern (SW) parts of M31 
show different kinematical properties, both sug- 
gest a rise in the outer part of the rotation curve 
(see Figure 5 and 6 in C09 and CIO, respectively). 

In this work, we show that magnetic dynamic 
effects constitute a clear and simple basis to 
interpret this feature. Some particular mod- 
els of magneti cally driven rotation curves have 
been presented (|Nelsonlll988UBattaner et aT]|l992t 
Battaner fc Floridoll995ll2000[l2007tlKutschera fc Jaloch 



two physical components: the dark matter halo 
and the magnetic field. 

2.1. Disc model 

T he stellar disc is assumed to be exponen- 
tial ( Freeman 1970l) . being the surface mass den- 
sity 

E exp(^], (1) 



RdJ 

where Eo is the central density of stars and Rd, 
the radial scale factor. 

The contributi on to the circular velocity of the 
stellar disc is (see Bin nev fc Tremaind 1987 ) 



v d (r) = 4TrGX R d y 2 [I (y)K (y) - h^K^y)], 



200-4 iTsiklauril 120081) . Although those models 
were originally developed to explain flat rotation 
curves without dark matter, our purpose here is 
more conventional and we will consider the con- 
tribution of both a NFW dark matter halo and a 
magnetic field added to match the velocity in the 
outermost region. 

Magnetic fields ar e known to slow l y decrease 
with radius (see e.g. iBeck et al. 19961 : Han et all 
19981 : lBeckll200lL 12004 l2005l )~and therefore they 
become increasingly important at the rim. For 
large radii, an asymptotic 1/r-profile for the field 
strength provides an asymptotic vanishing mag- 
netic force. This 1/r-profile has been found to 
match the polarized synchrotron emis sion of the 
Milky Wav (iRuiz-Granados et al.ll20l"oh and NGC 
6946 (|Beckl I2OO7I) . and will be considered in this 
work. 

2. Mass models and the magnetic field of 
M31 

In this section, we briefly present the luminous 
and dark mass models used to describe the ob- 
servational rotation curve of M31. As we are in- 
terested in the outer region (i.e. distances higher 
than 10 kpc), we follow C10 and we neglect the 
bulge contribution. We will also adopt from C09 
and C10 the parameters describing stellar disc and 
gaseous distribution, respectively. Therefore, we 
will only leave as free parameters those describing 



(2) 

where y = r/2i?^ and I , K , K\ are the modi- 
fied Bcsscl's functions for first and second order. 

The gaseous disc mainly contains HI, molecular 
gas and helium. The estimated gaseous distribu- 
tion of M31 represents approximately a 9% (for 
C09) and 10% (for C10) of the total mass of the 
stellar disc (see Table [I]). 

For the spatial distribution of the gas, we again 
assume the same exponential law given by Equa- 
tion ([J). 

2.2. Halo model 

The dark matter halo is descr ibed here by a 
NFW profile (jNavarro et al.lll996l ) 



Ph(r) 



ScPc 



1 



(3) 



where S c is a characteristic density contrast, 
p c = 3Hq/8itG is the critical density and Rh 
is the radial scale factor. The contribution to 
the circular velocity du e to this profile is given 
by (jNavarro et al.lll997l) 



l+cx 



1 ln(l + cx) 
V h (x) = V 200 \ , r 

V x M 1 + c ) - — 



(4) 



where V200 is the velocity at the virial radius -R200 
which is assumed ~ 160 kpc for C09 and -R200 ~ 
200 kpc for C10, c is the concentration parameter 
of the halo which is defined as c = -R200 / Rh and 
x is r/i?200- According to C09 and C10, the total 
dark halo mass at the respective virial radius is 
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-M200 ~ 10 12 M®. Our free parameters are V200 
and c. Values from C09 are c = 20.1 ± 2.0 and 
V 2 ao = 146.2 ± 3.9 km/s. Values from CIO are 
c = 12 and V 2 oo ~ 140 km/s. 

2.3. The regular magnetic field of M31 

The first measureme nts of the m agnetic field of 
M31 were obtained by iBeckl (jl982n using the po- 
larized synchrotron emission at 2700 MHz, show- 
ing that a magnetic pattern was aligned with HI 
structures and formed a ring at r ~ 10 kpc. 
They obtained a s t rengt h of B reg = 4.5 ± 1.2 /iG. 
Berkhuiisen et al.l ([20031 ) used observations of po- 
larized emission to show th a t M3 1 hosts an ax- 
isymmetric field. Han et al. ( 19981 ) used Faraday 
rotation measurements to show that the regular 
field could extend at least up to galactocentric 
distances of r ~ 25 kpc without significant de- 
crease of the strength and at least with a z ~ 1 
kpc of h eight above the g alacti c plane. However, 
recently, IStepanov et al. (2008) have shown that 
these results contained a significant contribution 
of Galactic foregrounds and so, it is difficult to in- 
fer the model from these measurements. In any 
case, as we show below, the detailed structure of 
the field is not relevant for the rotation curve. 

iFletcher et al. ( 2004 ) presented a detailed 
study of the regular field of M31 based on multi- 
wavelength polarized radio observations in the re- 
gion between 6 and 14 kpc. By fitting the observed 
azimuthal distribution of polarization angles, they 
found that the regular magnetic field follows an 
axisymmetric pattern in the radial range from 8 to 
14 kpc. The pitch angle decreases with the radial 
distance, being p ~ —16° for distances r < 8 kpc, 
and p ~ —7° for r < 14 kpc. This fact implies 
that the field becomes more tightly wound with 
increasing galactocentric distance. They found 
a total field strength (i.e. regular and turbulent 
components) of » 5 /xG. For the regular field they 
found that it became slowly lower, reaching at 
r ~ 14 kpc a strength of ~ 4.6 /xG. 

In this paper, our basic assumption is that the 
regular magnetic field of M31 is described by an 
axisymmetric model, that extends up to 38 kpc. 
For this model, the components in cylindrical co- 
ordinates are given by 



where p is the pitch angle and B$(r) is the field 
strength as a function of the radial distance. As 
shown below, from the point of view of the descrip- 
tion of the rotation curve, the relevant component 
is Brj,. Here, we will assume that p = 0°, as we 
are mostly interested in the outer region where p 
is very low. Several previous probes also indicates 
low values of p. In any case, this is not a strong 
assumption in the sense that a non-zero p value 
can be absorbed into the field strength (Bi) as a 
different amplitude. For the field strength, as a 
baseline computation, we shall consider a radial 
dependence of Bq(t) or equivalently B<p given by 



B l 



(7) 



where r± represents the characteristic scale at 
which B (r) decreases to half its value at the galac- 
tic centre and B\ is an amplitude in which we 
are absorbing the cos(p) factor. This expression 
has an appropriate asymptotic behaviour, in the 
sense that we obtain a finite value when r is close 
to the galactic center (r — > 0), and asymptoti- 
cally tends to ex 1 I r when r — ¥ 00, as suggested 
iBattaner fc Floridol (12007 ). Observations carried 
out by IFletcher et al.l (|2004[) established that 



B^r = Ukpc) 4.6/xG. 



(8) 



B r =B (r)sin(p) 
B^ = B (r)cos(p), 



(5) 
(6) 



This observational value at this radius was con- 
sidered as fixed in our baseline computation. By 
substituting into Equation ([?]) ■ we can find a rela- 
tion between B\ and r\ that allows us to re-write 
Equation ([7]) in terms of a single free-parameter 

(n) 

4.6n + 64.4 

B A r ) = ' ( 9 ) 

ri + r 

where r is given in kpc and B<p(r) in fiG. 
3. Methodology 

3.1. Observational rotation curve of M31 

We have considered the two datasets from C09 
and C10. They consist on a set of 98 and 29 mea- 
surements of the circular velocity (and their as- 
sociated error bars) respectively, which were ob- 
tained with the high-resolution observations per- 
formed with the Synthesis Telescope and the 26- 
m antenna at the Dominion Radio Astrophys- 
ical Observatory (C09) and with the wide-field 
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and high-resolution HI mosaic survey done with 
the help of the Westerbork Synthesis Radiotele- 
scope and the Robert C. Byrd G reen Bank Tele- 
scope (GBT) (jBraun et al.l l2009h ■ For our pur- 
poses, we consider only distances higher than 
r > 10 kpc to illustrate the effects of the mag- 
netic field. The actual data points on the rota- 
tion curve were obtained after fitting a tilted ring 
model to the data, and assuming a distance to 
M31 of 785 kpc. C09 derived a value of the incli- 
nation angle of i ~ 74°, which is lower than that 
derived from optical surface photome try measure- 
ments (|Walterbos fc Kennicuttlll987l ) and by CIO 
(i ~ 77°). The data points are plotted in Figure Q] 
in two separate panels. 

3.2. Influence of the magnetic field on the 
gas distribution 

The presence of a regul ar magnetic field af- 



fects the gas distribution (IBattaner et al.l 11992 ; 



Battaner fe Floridol [l99l l2000|). The jhiid mo 
tion equation can be written as (see e.g. IBattaner 
1996|) 



P^+MrWo+VP = nF+^B-VB-V ( f- 

Ot 47T \ 57T 

(10) 

where p is the gas density; v , the velocity of the 
fluid; P, the pressure; F, the total force due to 
gravity and B, the magnetic field. We assume 
standard MHD conditions i.e. infinite conduc- 
tivity. Equation (fTO]) is simplified by assuming 
axisymmetry and assuming pure rotation, where 
vb — (vor,vo<p,voz) = (0,9,0) even if these condi- 
tions are not necessary regarding the dynamic ef- 
fects in the radial direction. Taking into account 
all these facts, the motion equation in the radial 
cylindrical coordinate is 



where $(r) is the gravitational potential; F™ as , 
the radial component of the magnetic force, and 
P the pressure of the fluid. We can assume that 
pressur e gradients in the radial d irection are neg- 
ligible (Batt aner fc Floridol l2000h . In this case, 
then the radial component of the magnetic force 
is given by 




prnag = I _V , 



(12) 



and the contribution of the magnetic field to the 
circular velocity is given by 



B% , UB 2 



mag Airp \ r 2 dr ' 
3.3. Modelling the rotation curve 



(13) 



The rotation curve is obtained, as usual, by 
quadratic summation of the different contributions 

9{rf = v b {r) 2 +v d (r) 2 +v h (r) 2 +v rnag (r) 2 , (14) 

where we explicitly set v b (r) = as mentioned 
above. 

3.4. Model selection 

For the luminous mass models, the different pa- 
rameters are considered as fixed values in our anal- 
ysis (see Table [J). For the NFW dark halo, we 
constrain the V200 and c parameters, allowing a 
range for V200 between 100 and 220 km/s with 
steps of 0.5 km/s and c g [5, 30] with steps of 0.3. 
The contribution of the magnetic field to the ro- 
tation curve is fitted through one free parameter, 
r\ , that we are considering which is equivalent to 
fit £>0 as we discussed above. For this parameter, 
we have explored values in the range from 1 to 
1000 kpc. Our analysis is based on a reduced-^ 2 
as the goodness-of-fit statistic. Thus, the best-fit 
parameters are obtained by minimizing this func- 
tion 
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qvnodel \2 



X" 



N-M ^ 

i=l 



(15) 



where N is the total number of points to which 
we have measured the rotational velocity and de- 
pends on the considered dataset (N = 74 for C09 
and N — 27 for C10) and M is the number of free 
parameters. The sum runs over the observational 
data points, being 6° bs the observed velocity and 
Q-modei tne moc [ c u cc i velocity, which depends on 
the particular model. We shall consider two mod- 
els: one without magnetic contribution (DM) and 
another with the magnetic field (DM+MAG). Fi- 
nally, Ui is the observational error bar associated 
to each data point. 

4. Results and discussion 

Our main results are summarized in Figure [T] 
The dotted line shows our best-fit rotation curve 
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Dataset 


Disc parameters 


C09 


Rd = 5.6 kpc 
M d = 7.1 x 10 10 M Q 
M gas ~ 6.6 x 10 9 M Q 


CIO 


R d = 4.5 kpc 
M d = 8.0 x 10 10 Mq 
M 3QS ~ 7.7 x 10 9 M 



Table 1: Fixed parameters for bulge and disc. 



for C09 and CIO when considering only the usual 
dynamical components (the stellar component and 
the dark halo at this range of distances), while the 
solid line shows the result when adding up also the 
magnetic contribution. 

The most important resu lt is that, as shown 



Battaner fc Floridol (|2000h . the effects of mag- 



netism on the rotation curve are only relevant at 
large radii (in this case of M31, at distances larger 
than about 25 kpc). In this sense, the radial range 
of the datasets of C09 and CIO are optimal to ob- 
serve the magnetic effects. 

Table [2] summarizes our results for the best- 
fit solutions without (labelled as DM model) 
and with magnetic field influence (DM + MAG 
model) . As shown, magnetic effects on the gaseous 
disc significantly decrease the value of the reduced 
X 2 statistic for both datasets. Specially for C09, 
the fit is significantly improved when taking into 
account a new parameter (A% 2 = 6.5). The radial 
scale factor of the magnetic field (r*i) is uncon- 
strained in both cases, but shows clear preference 
for high values, which means that the best-fit so- 
lution for the field slowly decreases with the radial 
distance in the considered interval (i.e. between 
10 and 38 kpc). For example, at r ~ 38 kpc, the 
field is found to be > 4.4 fxG for C09 and 
B</> > 4.0 xtG for CIO for this best-fit solution. 
Both values are comp atible with the streng th of 
the field obtained by iFletcher et al. (2004) who 
found a nearly constant strength of the regular 
field of about ~ 5ixG between 6 kpc and 14 kpc. 
Moreover, when no radial variation of the strength 
is considered, (i.e. if r\ — > oo), and we fit for the 
amplitude B^, we obtain that B^ = 4.7±°,^ /zG 
which it is again compatible with results discussed 
above. This suggests that the data do not require 
an important radial variation of the strength of the 



field for the considered range of distances, or in 
other words, the contribution of the second term 
in the r.h.s. of Equation (fT2|) is negligible (i.e. 
dB^/dr <C B^/r). Therefore, if we had consid- 
ered another radial profile (e.g. exponential), we 
would have found a large radial scale factor too. 
The azimuthal component of the field is practi- 
cally constant between 10 and 38 kpc. 

Our results imply large magnetic fields at large 
radii. How these are produced lies beyond our 
scope. On the other hand, we would expect the 
extragalactic field to be also of t his order of mag- 
nitud e. Theoretical predictions ()Dar fc de Ruiula 



2005) suggest values of the level of few /zG for the 



intergalactic m a gnetic fields (hereafter IGMF). 
KronbersJ((l99i ); lGovoni fc Ferettil ( |2004h ; lKronberg| 



([20051 ) have reviewed observations of /iG level in 
rich clusters, though no direct measurements have 
been reported for the IGMF in the Local Group 
near M31. The observational evidence for IGMF 
is still weak, but quite strong IGMF near galaxies 
cannot be disregarded. 

We finally note the apparent discrepancy be- 
tween our derived parameters for the DM model 
and those obtained by C09 and C10 (see Sect. 12. 2} . 
However, it is important to stress that we are re- 
stricting the fit to the outer region of M31 (r > 
10 kpc). In this outer range, there is a weaker de- 
pendence of the shape of the rotation curve on the 
concentration parameter c, and thus lower values 
of c are found because the fit tries to compensate 
the rising behaviour in the outer part. The inclu- 
sion of the magnetic field corrects this apparent 
discrepancy, and in this case the values of c are 
now compatible with those obtained by C09 and 
C10. 
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Conclusions 



The rotation curv e of M31 (jChemin et al.ll2009t 
Corbelli et aDl2O10T) is rather representative of the 
standard rotation curves of spirals for r < 30 kpc, 
but it seems to rise out to, at least 38 kpc, the 
limiting distance of the observations. Indeed, this 
behaviour is not restricted to M31, and we are 
probably dealing with a common dynamical fea- 
ture of many other spirals. Therefore, the out- 
ermost rising rotation curve is a very important 
theoretical challenge. 

It is certainly a challenge as the standard dark 
matter halo models, in particular the universal 
NFW profiles, do not account for this dynamical 
unexpected behaviour. 

A conventional galactic model, with bulge, disc 
and dark matter halo, has been shown to provide 
good fit to the data in the range r < 20 kpc (C09, 
CIO). Here, we take advantage of these results 
and we do not fit any of the luminous components, 
which are taken to be exactly the same as those 
proposed by C09 and CIO. We have restricted our 
study to the region r > 10 kpc, and we only fitted 
the parameters describing the dark matter halo, 
and the magnetic field contribution. Our main 
conclusion is that magnetic fields are not ignor- 
able for explaining large-scale dynamic phenom- 
ena in M31, producing a significant improvement 
of the fit of the rotation curve at large distances. 
Moreover, the required field strength of the reg- 
ular component (B^ ~ 4 /iG) is fully consistent 
with the measured magnetic field in M31 at least 
up to r ~ 15 kpc. 

This conclusion seems very reasonable, as mag- 



Model 


Parameters 


C09 


C10 


DM 


V200 


(km/s) 


160.2 ±2.0 


132.1+g-J 






c 


12.3 ±0.6 


19-l±!i 






x 2 


19.8 


1.1 


DM + MAG 


V200 


(km/s) 


133.8^ 


120.01^ 






c 


22.7±i"f 


25.0+2'g 




r\ 


(kpc) 


> 888.0 


> 185.0 






x 2 


13.3 


0.6 



Table 2: Best-fit for the rotation curve with and 
without the contribution of magnetic fields for r > 
10 kpc. 



netic fields are amplified and act "in situ", and 
therefore they become increasingly important at 
the rim, where gravity becomes weaker. 

The best-fit model of the magnetic fields re- 
quires a field strength slowly decreasing with 
radius. This slow decrease is compatible with 
present values of the strength derived from ob- 
servations of the polarized synchrotron emission 
of the disc, but clearly we need measurements of 
Faraday rotation of extragalactic sources at this 
large radii to confirm that the magnetic field is 
present up to this distance and to trace unam- 
biguously the regular component. Hence, f uture 
exper iments such as LOFAR0 and SKA0 (Beck 
2009), will be extremely important, allowing a de- 



tailed explorations on the galactic edge as well as 
in the intergalactic medium. 
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Fig. 1. — Best-fit solutions for the rotation curve 
of M31, with and without including the contri- 
bution of a regular magnetic field component. 
Top shows the C09 dataset and bottom the CIO 
dataset. Asterisks and rombhus represent the ob- 
servational data with the associated error bars. 
The solid line is the best fit derived in this paper, 
including the contribution of the regular magnetic 
field over the gaseous disc. The dotted line is the 
best-fit model obtained without the contribution 
of the magnetic field. 
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